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ABSTRACT 

Context. The classical method to infer microturbulent velocity in stellar spectra requires that the abundances of the iron lines are not 
correlated with the observed equivalent widths. An alternative method, requiring the use of the expected line strength, is often used to 
by-pass the risk of spurious slopes due to the correlation between the errors in abundance and equivalent width. 
Aims. To compare the two methods and identify pros and cons and applicability to the typical practical cases. 

Methods. I performed a test with a grid of synthetic spectra, including instrumental broadening and Poissonian noise. For all these 
spectra, microturbulent velocity has been derived by using the two approaches and compared with the original value with which the 
synthetic spectra have been generated. 

Results. The two methods provide similar results for spectra with SNR>70, while for lower SNR both approaches underestimate 
the true microturbulent velocity, depending of the SNR and the possible selection of the lines based on the equivalent width errors. 
Basically, the values inferred by using the observed equivalent widths better agree with those of the synthetic spectra. In fact, the 
method based on the expected line strength is not totally free from a bias that can heavily affect the determination of microturbulent 
velocity. Finally, I recommend to use the classical approach (based on the observed equivalent widths) to infer this parameter. In 
cases of full spectroscopical determination of all the atmospherical parameters, the difference between the two approaches is reduced, 
leading to an absolute difference in the derived iron abundances of less than 0.1 dex. 
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' 1. Introduction 

One of the major limitations in the model atmospheres based 
. on 1 -dimensional geometry is an unrealistic treatment of surface 
' convection. Usually, the convective energy trans fer is computed 
I by a dopting the classical mixing-length theory jBohm-Vitensel 

■ 11958 ). a parameterization aimed to by-pass the lack of a rigor- 
' ous mathematical formulation for the convective energy transfer 
, On the other hand, this parameterization cannot describe realis- 
• tically the convective phenomena, ignoring their time-dependent 

■ and non-local nature. 

Basically, photospheric turbulent motions are divided in two ty- 
pologies, according to the size of the turbulent element: if this is 
[ smaller than the mean free path of the photons we call it micro- 
turbulence, otherwise macroturbulence. 

Such motions broaden the spectral lines, and the necessity to 
parametrize these effects arises from the discrepancy between 
empirical and theoretical Curves of Growth (COG) on their flat 
part. Microturbulent velocity (hereafter v,) is a free-parameter 
introduced to remove this discrepancy and to compensate for the 
missing broadening. Usually, v, is included during the spectral 
synthesis in the Doppler shift calculation, adding in quadrature 
Vt and thermal Doppler line broadening. 

Some general objections can be advanced on this simplification: 
( 1 ) the distinction between micro and macroturbulence is only an 
arbitrary boundary, because based on the uncorrected assump- 
tion that the velocity fields work on a discrete scale, while these 
velocity fields rather run on a continuous scale; 



(2) the 1-D model atmospheres (usually adopted in the chemi- 
cal analysis) include depth-independent microturbulent velocity. 
Observational hints point toward a different scenario, see e.g. the 
works about Arcturus bv lGravl (1198 1 ) and Takedal (119921) : 

(3) the turbulent velocity distribution is generally modeled as 
isotropic and Gaussian; this is contradicted by the observa- 
tion of the granulation on the solar surface; as pointed out 
by Nordlund et al. (199^) and |Asplund et al. (2000), the phe- 
nomenon of photospheric granulation is characterized by lam- 
inar (and not turbulent) velocities, due to strong density stratifi- 
cation; 

(4) mixing-length theory describes the general temperature 
structure of the photosphere, but it is unable to model some 
relevant convective effects, as the photospheric velocity field 
and the horizontal fluctua tions in the temperature structure 
dSteffen & Holwe^eiil2002h . 

The classical and empirical way to estimate v, is to erase any 
slope between the abundance A(Fe)[3 and the reduced equivalent 
width log[o(£'W//l) (hereafter EWR, where EW is the equivalent 
width), because v, works preferentially on the moderated/strong 
lines, while the weak ones are more sensitive to abundance. Even 
if the least square fitting of these two quantities is usually used 
to infer v,, this method requires that the two quantities have un- 
correlated errors. This condition is not respected in this case, 
because the errors in EWR and A(Fe) are obviously correlated, 



' I adopted the canonical formalism A(Fe) = \ogjg(N(Fe) /N(H)) + 
12. 
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with the risk to introduce a spurious slope in the A(Fe)-EWR 
plane. 

To byp ass this obstacle and provide a more robust diagnostic 
to infer Vt, iMagainI d 19841) proposed to derive v, by using the 
theoretical EW (hereafter EWT) instead of the observed one. In 
this paper, I investigated the impact of these two approaches to 
derive v,, in order to evaluate their reliability and to estimate the 
best one to use in the typical chemical abundance analysis. 



2. Artificial spectra: general approach 

The basic idea is to generate sets of synthetic spectra with a given 
and a priori known v, (with inclusion of instrumental broaden- 
ing and noise) and to analyze them with the two methods. Such 
a procedure allows to estimate the absolute goodness of the two 
methods (comparing the derived v, with the original one) but 
also the relative accuracy between the two approaches, in order 
to highlight systematic effects, if any. 

Different sets of artificial spectra have been generated with the 
following method: 

(1) synthetic spectra for a giant star with retf=4500 K, log g=2.0, 
[M/H]= -1.0 dex and v,-2 km/s have been computed with the 
spectro-synthesis code SYNTHE by R. L. Kurucz in its Linux 
version (Sbordone et al., 2004, 2005), starting from LTE, plane- 
parallel ATLAS 9 model atmospheres and by adopting the line 
compilation by R. L. Kurucz 0; the synthetic spectra cover a 
wavelength range between 5500 and 8000 A; 

(2) the spectra have been convolved with a Gaussian profile in 
order to mimic the instrumental resolution of the modern echelle 
spectrographs. A spectral resolution of Aid A- 45000 has been 
adopted (analogous to that of several spectrographs as UVES, 
CRIRES, SARG); 

(3) the broadened spectra have been re-sampled to a constant 
wavelength step of 5x- 0.02 A/pixel (similar to that of UVES 
echelle spectra); 

(4) poissonian noise has been added to the rebinned spectra, by 
using the IRAF task MKNOISE, in order to reproduce different 
noise conditions. SNR of 25, 40, 50, 70, 100 and 200 have been 
simulated. Due to the random nature of the noise, I resorted to 
MonteCarlo simulations, and for each value of SNR, 200 artifi- 
cial spectra have been generated and analyzed independently. 

I selected from the Kurucz linelist ~250 Fe I transitions in the 
spectrum, taking into account only unblended lines or those with 
negligible blending, similarly to what is done on real observed 
spectra. For all spectra, EWs have been measured with the auto- 
matic code DAOSPEC (IStetson & Pancinoll2008l) . 



3. Determination of v, and caveats 

For all the synthetic spectra described above, least square fits 
have been computed both in the A(Fe)-EWR and A(Fe)-EWT 
planes, choosing the value of v, that minimizes the slope. All 
the other parameters (Teff, log g and [M/H]) were kept fixed, 
according to the original ones of the synthetic spectra0. 



-4.5 



^ http://kurucz.harvard.edu/linelists.html 

^ Note that this is a simplifying assumption. In fact, the most adopted 
method to perform a chemical analysis (the so-called fuJJ spectroscopic 
method) is based on the following requirements: (i) the iron abun- 
dance is independent of (constraint for Teff); (ii) the same abun- 
dance should be derived from single and ionized iron lines (constraint 
for log g); (iii) the iron abundance is independent of the line strength 
(constraint for v,); (iv) the model metallicity well reproduces the iron 
abundance (constraint for [M/H]). 
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Fig. 1. Behavior of the true equivalent widths measured in the 
synthetic spectrum (without Poissonian noise) computed with 
the atmospheric parameters of a giant star, as a function of the 
quantity logn/g/) - O^ex- 



The most adopted proxy for the theoretical line strengthcan 
be written as logio(§/) - Oxex (see also Eq. 14.4 in lGravLTl992h . 
where logio(g/) is the oscillator strength, 9 the Boltzmann factor 
(defined as 6- 5040/reff) and Xex the excitation potential. Note 
that EWT does not include terms strictly inferred by the observa- 
tions but only related to the adopted model atmospheres (6) and 
to atomic data of the involved transitions (logi()(g/) andxex)- 

Fig. [T] shows the comparison between the true equivalent 
width computed for the synthetic spectrum of the giant star (be- 
fore adding noise) and the proxy used for EWT. As appreciable 
in figure, the quantity logio(g/) - Oxex is proportional to the true 
equivalent width of the synthetic spectrum (at least in the EW 
range shown in Fig. [1]), being a monotonic function of the true 
EW. In this work, I adopt the quantity logio(g/) - Oxex as the 
proxy for EWT because this is the formalism usually adopted 
in the chemic al analysis that follow the method described by 
IMagainI (II984I) . 

A threshold for the strongest lines has been applied: only 
lines weaker than log(EW//l)= ^.85 have been included in the 
analysis, coiTesponding to ~78 mA at 5500 A and ~II3 mA at 
8000 A. Stronger lines (basically located along the flat/damped 
part of the COG) have been excluded in order to not introduce 
additional effects (both observational and theoretical) that could 
affect the determination of v,. In fact: (i) strong features devi- 
ate from the Gaussian approximation (due to the development 
of Lorenztian wings), with the risk to under-estimate the EW 
(as discussed in Stetson & Pancino, 2008, see their Fig. 7), and 
(ii) these lines are affected by other effects, in addition to the mi- 
croturbulent velocity, i.e. hyperfine/isotopic splitting, damping 
treatment and modeling of outermost atmospheric layers: the in- 
clusion of such lines could lead to a different v,, very sensitive to 
the treatment of these effects for the iron lines, but it is not sure 
that the derived v, will be a ppropriate also for other elements 
(as discussed bv lRvanL Il998i) . Note that, in principle, stronger 
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lines can be taken into account in the abundance analysis (and 
for some elements like Ca or Ba only strong/damped transitions 
are often available), but they should be excluded from the deter- 
mination of V,. 

It is noteworthy that the approach used in this work repre- 
sents only a simple case with respect to the analysis of real spec- 
tra, because it is based on some simplifications: 

(i) large number of unblended lines. In real spectra, the num- 
ber of lines is limited by the quality of the available oscillator 
strengths and by the wavelength coverage; 

(ii) only v, is adopted as a free parameter to determine, because 
the other atmospherical parameters are known a priori from the 
computation of the synthetic spectra; 

(iii) the used spectra are flat and their continuum (virtually) set 
to the unity. Hence, distortions of the continuum due to residual 
fringing, uncorrected merging of the echelle orders or residual 
of blaze function are not included; 

(iv) this analysis is limited to features that can be well repro- 
duced by a Gaussian function, while in real cases some features 
with prominent damping wings can be included in the analysis 
by adopting a Voigt profile fitting (but taking into account the 
caveats discussed above for strong lines). 



4. Impact of SNR 

For each SNR, 200 artificial spectra were generated and ana- 
lyzed independently, and the mean value of v, of each group of 
SNR was derived. An important step in this kind of analysis is a 
careful selection of the lines used to infer v,, excluding features 
not well measured that could introduce a bias in the v, determi- 
nation. In order to quantify the impact of inaccurate EWs in the 
Vi estimate, I determined v, by considering different error thresh- 
olds and by taking into account only lines with errors ctew less 
than 10, 20, 30 and 50%; these errors have been computed by 
the DAOSPEC code from the standard deviation of the local flux 
residuals and represent a 68% confidence interval on the derived 
value of the EW (see Stetson & Pancino, 2008, for details). The 
discard of lines with inaccurate EWs have been performed indi- 
vidually for each artificial spectrum. As sanity check, the EWs 
measured in a synthetic spectrum with SNR=25 have been com- 
pared with the true EWs of the original spectrum (before adding 
noise). The results are shown in Fig.|2j the inner panels show the 
behavior of the relative error in percentage as a function of the 
true and measured EW (left and right inset, respectively). These 
plots confirm the reliability of the measured EWs and show that 
the errors in the EW measurement affect mainly the weak lines, 
while the strong (>40 mA) lines are well measured with small 
uncertainties. 

The results concerning the determination of v, with the two 
methods for the sample of synthetic spectra described above are 
plotted in Figl3] where empty and filled circles indicate the val- 
ues of V, derived by using the EWR method and EWT one, re- 
spectively. The errorbars indicate the dispersion around the mean 
of each set of simulations of a given SNR. For sake of clarity, the 
original value of v, of the synthetic spectra is plotted as a dashed 
line. Typically, the dispersion around the mean in each group 
ranges from a few hundredths of km/s (for SNR= 100 and 200) 
up to 0.2 km/s for SNR= 25 and considering lines with resid- 
ual errors less than 50%; note that such dispersions are basically 
compatible with the corresponding errors of the computed slopes 
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Fig. 2. Difference between the measured and true EW in a syn- 
thetic spectrum with SNR=25 as a function of the true EW. In 
the left inset is shown the behavior of the error in percentage of 
the measured EW as a function of the true EW, while in the right 
inset is shown the same behavior as a function of the measured 
EW. 
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* For a sub-sample of synthetic spectra, the analysis has been re- 
peated by using the true EWs instead of EWT as defined previously. 



Fig. 3. Behavior of the measured microturbulent velocity as a 
function of SNR, derived wi th the classical method (empty cir- 
cles) and the prescription bv lMagai"nl d 19841) (filled circles) for a 
giant star. Each point is the average value obtained by using 200 
synthetic spectra of a given SNR and the errorbars indicate the 
dispersion around the mean. Dashed lines indicate the original 
value of V, of the synthetic spectra. Each panel shows the results 
obtained by considering different thresholds for the residual er- 
rors in the EW measurement. 



4 



A. Mucciarelli: Microturbulent velocity in stellar spectra (RN) 



In case of relatively high (>70) SNR, the two methods pro- 
vide very similar results and agree with the v, of the original syn- 
thetic spectrum. In fact, increasing the SNR, the errors in the EW 
measurements are reduced and the discrepancy between the two 
approaches is erased. On the other hand, spectra with SNR<70 
show V, always lower than the original value and the discrepancy 
increases when also lines with high uncertainties are taken into 
account. Generally v, inferred by using EWT turns out to be sys- 
tematically lower than that obtained with the classical method 
by 0.1 - 0.2 km/s, approximately. The reasons of this systematic 
under-estimate will be explained in Sect. |5] 
These findings highlight also the impact of measured EWs with 
high uncertainties: the inclusion of lines with EWs very uncer- 
tain (generally, the weakest lines in spectra with low SNR) wors- 
ens the discrepancy with the true v, of the original synthetic spec- 
tra. Thus, a careful selection of the used lines, considering the 
accuracy of their EWs, is crucial in the derivation of v„ in order 
to exclude additional bias. 

A simple analytical model can be derived by fitting the data 
shown in Fig. [3] The absolute value 5,,, of the difference be- 
tween the derived v, and the true one can be expressed as 
(5,,, - fa-Ew 1^ where fa-^w is a linear function of the adopted er- 
ror threshold ctew (expressed in %). The fit of the present dataset 
provides 6y^- (0.14 ■ ctew + 2.64)/S NR for v, derived by using 
EWR and 5,,= (0.19 ■ ctew + 5.16)/SNR when EWT is adopted. 
I repeated these checks also with synthetic spectra computed 
with the typical parameters of a dwarf star (reff=6000 K, 
log g=4.0, [M/H]= -1.0 dex and v,-l km/s), finding basically 
the same result. 
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Fig. 4. Behavior of the iron abundance as a function of EWR 
(upper panel) and EWT (lower panel) for a synthetic spectrum 
computed with v,-2 km/s. Filled points show the shift of a given 
line in the two cases due to the EW error. 



5. Is the EWT approach free from biases? 

An useful example to clarify the systematic offset between the 
two methods is to measure EWs of the red giant synthetic spec- 
trum before adding noise: all lines, regardless of their strengths, 
provide the same abundance. These EWs have been varied ran- 
domly by +lcr, where cr was computed with the classical for- 
mula by ICavrell (Il988h (cr = 1.5 ■ (SNR-^) ■ ^JFWHM ■ 6x, 
adopting SNR= 25, the wavelength step Sx- 0.02 A/pixel and 
a FWHM of A/R, coiTesponding to 0.12 A at 5500 A and 0.18 
at 8000 A) and the analysis repeated as described above. The 
error is assumed the same regardless of the size of the EW. 
Fig. |4] shows an example of this test, with the distribution of 
the Fe I lines in the A(Fe)-EWR and A(Fe)-EWT planes (upper 
and lower panel, respectively). Filled points show how a generic 
weak line (with initial A(Fe)= 6.5 dex) moves in the two dia- 
grams because of the effect of EW errors. In the upper panel of 
Fig. |4]the effect due to the correlation between the error in EW 
measures and derived abundances is clearly visible, while the 
correlation is avoided by definition in A(Fe)-EWT plane: this is 
the main reason behind the use of EWT. 

The mean slope, computed by taking into account 200 
MonteCarlo simulations of random error distribution, is of 
H-0. 13+0.02 in the A(Fe)-EWR plane, and -0.01±0.02 in the 
A(Fe)-EWT plane (ever keeping fixed the v, value at the orig- 
inal value of the synthetic spectrum, v,= 2 km/s). Such slopes 
point out the need to increase v, (of at least 0.2 - 0.3 km/s) 
to erase the residual slope when v, is optimized by using the 
EWR. On the other side, adopting the EWT, no relevant change 
in Vt are necessary. This finding reproduces exactly the bias dis- 



The results are very consistent with those obtained by using EWT, in 
agreement with the discussion of Sect.[3]and the finding of Fig.[T] 
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Fig. 5. Same of Fig. |4]but highlighting the lines detectable at 
SNR= 25 (grey points). Grey regions indicate the location of the 
undetectable lines. 



cussed bv lMaeailil d 19841 ) (and shown in his Fig. 2) but it seems 
to be contrary to the findings shown in Fig. [3] where the two 
methods provide lower v, with respect to the true value. In fact, 
Magain (1984) demonstrated that v, is systematically overesti- 
mated when derived with the classical method. 
How can we explain such a discrepancy? In the real case (and 
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also in the test performed with synthetic spectra, see Sect. |2]l, 
not all the lines can be measured and used to infer v,. The ef- 
fect of the slanting movement of the lines shown in Fig|4](upper 
panel) is critical for very weak lines (where the errors are more 
relevant) but at low SNR not all these very weak features are 
detectable, because they are too noisy or indiscernible from the 
noise envelope. Hence, many weak lines that contribute to the 
positive slope in the A(Fe)-EWR plane will go out of the de- 
tectability threshold, reducing the positive slope. 
For a given SNR, a reasonable and generally adopted threshold 
to measure spectral lines is to assume 3 times the uncertainty 
provided by the Cayrel formula. Thus, for SNR= 25 we can mea- 
sure absorption lines of at least 9 mA and 1 1 mA at 5500 and 
8000 A, respectively. Because the difference is negligible along 
the entire coverage of the spectrum, I adopted conservatively a 
threshold of 1 1 mA. Fig. |5]reports the same results of Fig. |4]but 
highlighting the lines really detectable assuming such a thresh- 
old (grey points), while the lines located in the grey shaded re- 
gions are excluded from the analysis. 

The mean slope of the 200 MonteCarlo simulations, com- 
puted on the surviving lines, turns out to be of -0.10+0.02 in the 
A(Fe)-EWR plane, and -0.14±0.02 in the A(Fe)-EWT plane. 
When we take into account the effective possibility to detect very 
weak lines with large uncertainties, the situation turns in to the 
opposite of the case described in FigH) 

- the threshold in the A(Fe)-EWT plane results to be diagonal, 
introducing a negative slope and the need to lower v, to erase 
any trend 0: 

- in both cases the derived v, will be lower than the original 
value of the synthetic spectrum (the same finding discussed 
above and shown in Fig. |3]l. 

In principle, all these effects are valid also for higher SNR, 
but become negligible because the quoted uncertainties in the 
EW measurement are of only a few percent. 

This experiment clarifies the origin of the systematic under- 
estimate of V, with the two approaches (shown in Fig.|3]l. When 
V, is derived from EWR, the correlation between the EWR and 
abundance errors moves the lines in an asymmetric way on the 
A(Fe)-EWR plane (upper panel of Fig. |4]l, but the threshold of 
detectability of the lines will populate mainly the region with 
A(Fe)>6.5 dex with respect to other one (as appreciable in the 
upper panel of Fig. |5]l, leading to a negative slope (and an under- 
estimate of V,). 

When V, is derived from EWT, the threshold of the line selec- 
tion shown in the lower panel of Fig.|5]increases the asymmetric 
distribution of the lines in the plane, ever leading a negative slope 
(as demonstrated by the average slopes in the two cases) and a 
slightly lower v, value. 

Obviously, the inclusion of weak lines with inaccurate mea- 
sures will worsen the discrepancy (lower panels in Fig.|3]l. 

Summarizing, v, inferred by using the A(Fe)-EWT plane are 
not totally free from observationaJ biases. In fact: 
(a) the use of A(Fe)-EWT plane partially erases the effect due to 
the correlation between errors in EW and abundances, but also 
introduces another bias due to the threshold of detectability of 
the lines (more relevant and crucial at low SNR), as described 
above; 

^ Note that the slope of the boundary between detectable and un- 
detectable lines is roughly orthogonal to the shift of a line due to the 
correlation between the errors in EW and abundance. Basically, the two 
slopes have opposite sign. 



(b) errors in A(Fe) and EWT are still not totally independent 
from each other, because both quantities suffer in uncertainties 
from the adopted logio(^/) values. In fact, an overestimate of 
the oscillator strength for a given line will provide an overesti- 
mate of EWT but also an underestimate of A(Fe), introducing a 
negative, spurious slope when EWT is used to infer v,. Note that 
this effect is secondary, if accurate laboratory oscillator strengths 
are taken into account. For instance, th e typical errors of the Fe I 
log[Q(^/) in the critical compilation bv lFuhr & W iese (2006) are 
less than 25%. Such errors translate in an iron abundance error 
of the same order of magnitude obtained with an error of EW 
of ~25% or less. The correlation between the errors in iog^Qigf) 
and abundances introduce a spurious slope in the plane A(Fe)- 
EWT, in a similar way to the spurious slope in the A(Fe)-EWR 
plane due to the correlation between abundance and EWR dis- 
cussed above. Otherwise, one needs to consider also this sec- 
ondary bias only if inaccurate/inhomogeneous atomic data for 
iron lines are included in the linelist; 

(c) in the case where the temperature has been inferred by impos- 
ing the excitational equilibrium (erasing any slope in A{Fe)-Xex 
plane), also will be affected (in a stronger way for low SNR 
spectra) by the errors in EW and therefore also EWT will be af- 
fected by the same bias, by definition. 

Note that these last two aspects have not been deeply investi- 
gated in the present work 

6. Impact of the number of lines 

Previous results are based on a large number of Fe I lines (250), 
allowing to well sample both weak and moderately strong fea- 
tures. Such an assumption allows to avoid spurious slopes (re- 
gardless to the adopted method) due to inadequate sampling of 
line strengths. 

Some mid/high-resolution multi-object spectrographs offer 
wavelength coverages that limit the number of available Fe I 
lines. 

A simple test has been performed in order to clarify how 
the number of lines is critical. 50 Fe I lines have been extracted 
randomly from the measured lines in a red giant synthetic spec- 
trum with noise SNR=25 and v, was redetermined with the two 
methods. Such a procedure has been repeated by using 200 
MonteCarlo extractions. The average values of the v, distribu- 
tions is very similar to those obtained above but with a larger 
dispersion (cr ~0.3-0.4 km/s) with respect to that obtained in 
the case described in Section. H] 

If the same check is performed with a synthetic spectrum with 
SNR=100, the final distributions share dispersions around the 
mean of about cr ~0. 15-0.2 km/s, in both cases, while in the 
case discussed in Section |4] the dispersion was very small (see 
Fig.O. Also in this case, such dispersions are larger than those 
inferred from EWR and EWT when a large number of transitions 
is used. 

The most discrepant v, values in the derived distributions are 
obtained when an inadequate sampling of the weak and strong 
lines is adopted. In fact, if the experiment is repeated imposing 
the condition that the 50 random lines are extracted by sampling 
equally the entire EW range (thus providing a statistically signif- 
icant sample of weak and strong lines), the dispersion around the 
mean of the derived distribution of v, is reduced by about a factor 
of 1.5 and 2 for SNR=25 and 100 respectively. In this case, the 
distributions show dispersions around the mean that are barely 
consistent with those obtained with the entire linelist. I can con- 
clude that the number of used lines and the relative sampling 
between weak and strong lines are crucial elements that can af- 
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feet the determination of more severely than the choice of the 
EWR or EWT method in the analysis. 



7. Impact of the continuum placement 

As akeady pointed out above, the test performed with synthetic 
spectra neglects uncertainties arising from the continuum loca- 
tion, because the continuum in artificial spectra is set at unity 
and totally flat, while in real spectra several additional eff'ects af- 
fect its correct identification. In order to clarify how much this 
source of uncertainty in the EW measurements impacts v,, I per- 
formed two simple tests by considering synthetic spectra with 
SNR= 25: 

(1) the EWs have been remeasured by adopting a higher and 
lower continuum level, with respect to that derived in the previ- 
ous experiment. The variations in the continuum level are of +5 
and ±10%, translating to a systematic (positive or negative de- 
pending by the sign of continuum variation) of fset of the derived 
EWs (the same finding has been discussed in S tetson & Pancind 
(12008 ). see their Fig. 2). The differences between v, derived with 
the EWR and EWT methods are similar to those obtained with 
the original continuum level; 

(2) the synthetic spectra have been multiplied with a sinusoidal 
curve, in order to mimic a fringing-like pattern, similar to that 
observed in real spectrum. Also in this case, no relevant dif- 
ference (within the uncertainties) with the original analysis is 
found. 

The continuum tracement (at least in case where the con- 
tinuum is defined globally along the entire spectrum) does not 
affect dramatically the results described in Sect.|4]about the rela- 
tive difference between v, derived with EWR and EWT methods. 



8. A real case: UVES spectra of LMC giants 

In order to show a case based on real spectra, I derived v, 
with the two methods for s ome Large Magellanic C loud giants 
stars observed with UVES ( iMucciarelli etali |2009|) . with typi- 
cal SNR= 5 and originally analyzed following the prescriptions 
bv lMa^ainl(fT984l) . 

Microturbulent velocities computed with the classical 
method (and fixing the other parameters) are higher b y ~0.2 - 
0.25 km/s than those obtained with the prescription by Magain' 
([T984) : this difference translates to an iron abundance ~0.15 dex 
lower when the classical approach is adopted. 

A point to recall here is that the Tetf and v, can be correlated if 
these quantities are derived with the classical full spectroscopic 
method. In fact, the majority of weak lines have high;i^ex, while 
several strong transitions have low excitation potentials, leading 
to a correlation between ;^^ex and the strength of the lines. Also, 
variations in Teff and change differently single and ionized iron 
lines abundances, needing to re-adjust the surface gravity. The 
higher v, computed with the classical method implies the need to 
adjust spectroscopically both Tgff and log g, by increasing them 
and hence reducing the iron abundance difference (to less than 
0.1 dex). 

Summarizing, the absolute difference in iron abundance 
computed witli the two approaches is reduced (less than 0. 1 dex) 
if a full spectroscopic determination of all atmospherical param- 
eters is performed. Therefore, if one or more parameters are 
fixed to a certain value, this compensation does not take place 
and the absolute difference can be higher. 



9. Determination of v^: the best strategy 

I propose an easy strategy to minimize the effect of the bias dis- 
cussed above, deriving the v, from EWR, being the approach that 
gives the better results with respect to the true v, of the synthetic 
spectra used in Sect. |4| 

(i) being the determination of v, strongly affected by the cor- 
related errors in abundances and EWs, the lines with high un- 
certainties in the EW measurement need to be discarded. In fact, 
as demonstrated in Fig. [3] the inclusion of lines with inaccurate 
EWs leads to highly uncorrected v, with respect to the true value. 
Basically, the use of lines with EW errors less than 10% provides 
a good agreement with the v, value of the synthetic spectra. The 
reader can use the relations provided in Sect. |4] to estimate the 
amplitude of the bias in v, as a function of the adopted (Tew 
cutoff. Note that the computation of the EW uncertainty is not 
provided by all the codes developed to measure EWs, but a stan- 
dard error about the goodness of the fit for each measured line 
is recommended, through the standard deviation of the local flux 
residuals after the Gaussian fitting (as performed by DAOSPEC) 
or by adopt ing other formulae, i .e. the Cayrel formula or that 
proposed bv ' Ramirez et al.l (|200ll) : 

(ii) as shown in Fig. [3] even if an adequate lines selection is 
applied, a residual discrepancy between the derived and true v, 
of ~0.15 km/s remains in the most critical case of SNR= 25. A 
possible method to reduce this residual discrepancy is to per- 
form a linear fit obtained taking into account the uncertain- 
ties in both EWs and abundances (following the approach by 
Pre ss et al.l (1 19921) as implemented in the FORTRAN routine 
FITEXY). Errors in EWs have been estimated as described in 
Sect. |4] while the error in the abundance of each individual line 
has been derived by computing the abundance with the EW var- 
ied of ilcrgvy. This method to perform the linear fit allows to ob- 
tain values of v? consistent within the errors with the true value of 
the synthetic spectra (see the first panel of Fig.|3]l, reducing the 
residual discrepancy. Usually in the chemical analysis, the slope 
in A(Fe) — EWR plane is computed neglecting the uncertainties 
in the two variables (because this procedure needs the computa- 
tion of the error in both EW and abundance for each measured 
line) but this approach (or other approaches to find the linear re- 
lationship in data with unc ertainties in both directions, see e.g. 
iHogg. Bow & Langll20Toh is recommended. 

10. Conclusions 

I draw the following results from the above test: 

(1) the method to infer the microturbulent velocity from the 
EWT provides the same results (within the quoted uncertainties) 
of the classical method in case of moderate/high SNR (>70); 

(2) in case of SNR<70 both methods underestimate the micro- 
turbulent velocity and when low (25) SNR spectra are analyzed a 
better agreement with the true v, of the synthetic spectra is found 
with the classical method based on EWR; 

(3) the discrepancy with the true v, value (regardless of the 
adopted method) worsens when EWs with high uncertainties 
(generally, the weakest lines in low SNR spectra) are taken into 
account. This finding points out the need to estimate the quality 
of the measured EWs for each individual line (a procedure usu- 
ally not performed in the chemical analysis), in order to reduce 
this systematic underestimate; 

(4) these results are valid for both giant and dwarf stars, under 
the assumptions that the strong lines that deviate significantly 
from the Gaussian approximation are excluded from the analy- 
sis; 
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(5) the number of used lines and the adequate sampling of weak 
and strong Unes in the case of a small number of available iron 
lines affect the determination of Vt more than the adopted ap- 
proach to infer this parameter; 

(6) the spectroscopic optimization of all the atmospheric param- 
eters simultaneously is recommended, in order to compensate 
errors in the determinations of the parameters. 
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